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1. Introduction 

One of the most pressing issues at current and future hadron coUider experiments 
is the discovery of one (or more) Higgs bosons. The Standard Model (SM) predicts 
a single Higgs boson whose couplings to fermions are proportional to the mass 
of the fermion, hence the production of the Higgs boson with light quarks (e.g. 
bottom quarks) is suppressed compared to other production channels involving top 
quarks or heavy gauge bosons. However, in some extensions of the SM (such as the 
Minimal Supersymmetric Standard Model, or MSSM), the Yukawa couplings for b 
quarks can become strongly enhanced. In this scenario, the production of a Higgs 
boson with bb can dominate over other production channels and could provide a 
unique opportunity to directly probe the b quark Yukawa coupling. 

The theoretical prediction of bbh production at hadron colliders involves several 
subtle issues and depends on the number of b quarks identified, or tagged, in the 
final state. Potentially large logarithms arise from the integration over the phase 
space of the final state b quarks. These large logarithms appear at all orders in per- 
turbation theory and could spoil the convergence of the perturbative expansion of 

*Talk given by C.B. Jackson 



2, 2008 1:0 WSPC/INSTRUCTION FILE DPF2004jackson'bbh 



2 Dawson, Jackson, Reina and Wackeroth 



Vs=2 TeV 
M„=120 GeV 
|ii,=m^+M^/2 



^NLO.OS 
^NLO.US 
^LO.OS 
^LO.MS 



>20 GeV 
|ri|<2 



Vs=14TeV 
Mj=120 GeV 
|i„=mj+My2 



"NLO.OS 
*^NLO,MS 



Pj >20 GeV 
|r|<2.5 



0.3 0.5 0.7 1 



Fig. 1. ctnlo cr£,Q for pp,pp — > feb/i with two high-py 6 jets at y/s = 2 TeV (left) and y/s - 
14 TeV (right) as a function of the renormahzation/factorization scale fi, for Mjj = 120 GeV. 



total and differential cross sections. Currently, there are two approaches to calculat- 
ing the theoretical prediction of bbh production. Working under certain kinematic 
approximations, and adopting the so-called five-flavor-number scheme (5FNS), the 
logarithms can be resummed by using a bottom quark Parton Distribution Func- 
tion^'EI. Alternatively, working with no kinematic approximations, and adopting 
the so-called four-flavor-number scheme (4FNS), one can compute the cross section 
for pp, pp bbh at fixed order in QCD without resumming higher order coUinear 
logaritmsEl ■ S. 

In this talk, we present the results of the total cross section calculation for 
pp,pp — > bbh at next-to-leading order (NLO) in QCD for two different scenarios: two 
high-pT b jets and one high-pr b jet corresponding to two or one b jet(s) identified 
in the final state. Results for the case of no high-pT b jets have been reviewed in a 
previous study^l. For the case of one high-p-r b jet in the final state, we compare our 
fixed-order result (4FNS) with the prediction obtained by using gb — > bh (5FNS), 
when NLO QCD corrections have been included in both cases'^. 



2. Results for bbh Production 

In order to select the events with high-p^ b jets, we place kinematic cuts on the 
transverse momentum, > 20 GeV, and pseudorapidity, \rif,i\ < 2 (2.5) at the 
Tevatron (LHC), of the final state b quarks. Results from the case of two high-px b 
jets are shown in Fig. 1^. The two sets of curves correspond to different renormal- 
ization schemes for the bottom quark mass in the bottom quark Yukawa coupling. 
In both cases, the dependence on the renormalization and factorization schemes 
is greatly reduced by the inclusion of NLO QCD corrections. It is interesting to 
note that the MS scheme exhibits better behavior at small scales and also seems 
to improve the convergence of the perturbative series. The calculation of the NLO 
cross section for bbh for the case of (at least) one high-px b jet has been performed 
in both the 4FNS ^ El and the 5FNS El^ where the leading order process becomes 
gb(b) b{b)h). In a previous studyEI, the two cross section calculations were found 
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Fig. 2. Total NLO cross section for pp,pp — > bbh with at least one high-pr b jet at the Teva- 
tron(left) and LHC(right). The uncertainty bands are obtained by varying the renormalization and 
factorization scales independently around the central value {fir = l^-f = A'o/2). The solid(dashed) 
curves correspond to the 4FNS{5FNS) calculations. 

to be in agreement (within error), but the 5FNS calculation seemed to predict a 
slightly larger value for most of the Higgs mass range. In a subsequent paper ^, we 
have investigated the inclusion of diagrams containing top quark loops, that were 
previously neglected in the 5FNS, by implementing them into MCFM"!. The contri- 
bution from the top loop diagram(s) to the total cross section in the SM calculation 
can become numerically important. In fact, inclusion of these diagrams was found 
to reduce the 5FNS cross section by 15%(10%) at the Tevatron(LHC) over most of 
the Higgs mass range Results are shown in Fig. 2 for the comparison between the 
4FNS and 5FNS schemes after the inclusion of the top loop diagram in the latter 
calculation. Clearly, within the theoretical uncertainties, the SM calculations in the 
4FNS and 5FNS are in very good agreement. The agreement is preserved in the 
MSSM with large tan /3, since the top loop diagrams become negligible. 
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